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1. Introduction 

It is currently believed that most of the galaxies host supermassive black holes (SMBH) 
in their cores (2D]- Even our own galaxy host a black hole of almost 4 million solar 
masses |2nj in its center. In this scenario, the most likely event is that stars and other 
compact objects in the neighborhood of the SMBH will inspiral-down and plunge toward 
this central massive black hole. 

From the Theoretical point of view the two body problem in General Relativity 
remains one of the more interesting problems still unsolved. While comparable mass 
binary black holes are the most efficient astrophysical system regarding emission of 
gravitational radiation, they represent also the most challenging problem for numerical 
relativity. Current estimates show that during the final plunge of two nearly equal mass 
black holes the system radiates as much energy as in the whole previous inspiral history 
of the binary d El 121 13] ■ 

For the Laser Interferometric Space Antenna (LISA), the small mass ratio regime 
is particularly interesting since numerous sources will be 'visible' with huge signal to 
noise ratios [T^]. This scenario is tractable using perturbation theory of black holes, 
but proved to be a more challenging problem than initially expected, until relatively 
recently, when decisive progress toward the complete solution to the problem has been 
made, some of which we will briefly report here. 

2. The problem and its solution 



The theory to compute radiation from an orbiting particle around a black hole has 
been developed in the 70's [23 I2H 121] ■ Soon after it was realized that correcting the 
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geodesic of the particle by its radiation formally lead to infinities. In fact, if we try a 
straightforward solution to the case of a particle represented by a Dirac delta in orbit 
around a black hole we can write the field equations of general relativity as 

G^[h aP } = T llv = m J u^u u 5[x a - Xp (r)]/ dr (1) 

where h a p represents the perturbations around the larger black hole background metric, 
gjiJ (with determinant g). Xp(r) is the (parametrized by the proper time of the 
particle, r) coordinate location of the particle (with mass m) orbiting the large hole 
and = dx^ / dr is its four- velocity. 

Once we solve for the perturbed metric we can correct the trajectory of the particle 
by requesting it to be a geodesic of the background plus the perturbed metric 

The problem arises when trying to evaluate the connection term T^p on the trajectory 
of the particle since it can be readily shown that the metric behaves as 

mH uv , , 

h/jtu ~ , (3) 



\j (9aP + UaUp) (x a - X^ij)) 



T. 



i.e., diverging as the inverse of the distance to the particle. So, the small perturbations 
assumption breaks down near the particle. 

In order for us to extract any physically meaningful result this problem needs to be 
regularized. This proved to be a very challenging problem for over 25 years until 1997 
when Mino, Sasaki and Tanaka [21] published their seminal work on the computation of 
the self-force a particle traveling in a curved spacetime exerts on itself. They computed 
this by using two alternative methods. First using matched asymptotic expansions of 
the perturbed background of the larger black hole by the presence of the particle and the 
spacetime around the small particle tidally distorted by the larger hole. By matching 
this two expansions in a buffer zone of radius d around the particle, with m « d « M, 
one finds that requiring the consistency of the two expansions leads to the equations of 
motion. The second derivation is done via a Brehme-DeWitt approach by deriving a 
'conserved' rank-two symmetric tensor and integrating its divergence over the interior 
of the world tube surrounding the particle to obtain the equations of motion. 

The self-force (defined as the correction to the background geodesic) takes then the 

form 

FU (*) = mV% = F^ dlTect + (4) 

where F^ irect is computed from contributions that propagate along the past light cone 
and Ff ail has the contributions from inside the light cone, product of the scattering of 
perturbations due to the motion of the particle in the curved spacetime created by the 
larger black hole. 
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The self- force is then computed by taking the limit F^ = Ff ail (x — > x p ) where 

Ku = {^W* ~ V p G- a % & , - \u^V x G- ap A U^U S 'dT>(5) 

and where G~p , 5 , is the retarded Green's function. Primed indices denote evaluation 
at the particle's location while unprimed denotes evalutation at the nearby field point 
x (u™^ are obtained by parallel propagation.) Finally, r_ represents the proper time 
of a point on the worldline that intersects the particle's past lightcone of the field point 
x. As a result, one finds that the divergent terms mentioned above do not affect the 
motion of the particle. 

Independently, Quinn and Wald [22] come up with the same equation of motion (jSJ) 
(usually referred to as 'MiSaTaQuWa' force), but using an axiomatic approach: 

i) comparison axiom: Identifying the four-velocity w M and the four-acceleration at 
points Pi in a manifold (Mi,g$) and P2 in a manifold (M2,g$) via Riemann normal 
coordinates one find that the difference between forces is given by 

%r%) = h { (l v wC - v (d/3 ^ - Qv^g - v m h$ ^ } ( 6) 

ii) Flat spacetime axiom: If (M, g^ v ) is the Minkowski spacetime the perturbed metric 
is the half advanced and half retarded solution and the force vanishes 

hap = \ [Kp + Kp] , and r = mu a V^ u» = 0. (7) 

It is worth mentioning here a newer approach ^3] to the problem which decomposes 
the full fields into a 'regular' (which is different from the 'tail' part) and a 'singular' 
part. It can be proved that only the regular part contributes to the self-force, producing 
the same 'MiSaTaQuWa' result, and that the regular (metric) fields satisfy the vacuum 
equations of General Relativity (with no source terms due to the presence of the particle). 

The self-force calculation thus provides a local correction to the trajectory of the 
particle taking into account radiative and nonradiative effects. For a global method that 
takes into account the radiative effects only see, for instance Ref. (TB], that makes use 
of the energy and momentum balance of the radiation emitted to infinity and onto the 
black hole horizon to correct the particle trajectory over orbital cycles. 

3. The problems of the solution 

The proof that the equation of motion 

/ d 2 x» (0)M dx a dx?\ _ 

m \-*FT + ~d^^)- F ^ (8) 

leads to the correct, finite description of the particle due to its own field represent an 
enormous leap forward in our comprehension of the effects of gravitational radiation. 
Yet, expression (J3J for the self-force is formal and it was not immediately ready for 
implementation. It took a couple of years to develop the formalism that allowed 
explicit computations. The key observation, arrived to simultaneously by several authors 
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[El HUE], was to tackle the nonrotating, i.e. Schwarzschild black hole, background case 
first to make use of its spherical symmetry to expand the force into spherical harmonics. 
For each multipole (£,tcl), the full force is finite, and the divergent nature of the problem 
only reveals upon summing over (£,m). 

In the mode sum regularization method the force can be decomposed as (wehere 
we performed alrady the sum over the modes m) 

F self = J™ E - F ^)] (9) 

which upon addition and subtraction of A^ L + + C M / L (which leads to divergences 
upon summing over £) 

F seif = E [twM -(A»L + B» + C/L) 

- E -(^l + b» + c/l)] (io) 

t 

where for convenience we defined L = i + 1/2. The terms in brackets now lead to 
well defined, convergent series. The sum over the last bracket is a remainder normally 
denoted as . 

In a notable piece of work, the regularization parameters A^,B^,C^ and _D M 
have been computed analytically, first for the Schwarzschild [7j and then for the Kerr 
backgrounds [TT 



The missing part is now the full force, Ff ull (x p ), that has to be computed by 
numerically integrating the (perturbed) General Relativity field equations to obtain 
Ff ull (x). For a nonrotating large black hole, we can just solve the wave equations for 
the perturbations traveling on the curved Schwarzschild background, i.e. the Zerilli and 
Regge- Wheeler equation for the two polarizations of the radiation. We then know how 
to reconstruct the metric in the Regge- Wheeler gauge jTTj. 

We successfully followed this procedure in the case of a particle starting from rest 
at a finite distance from the large black hole |2J. In addition to the exact (numerical) 
computation of the self-force, we obtained an analytic approximation 

f.w = ^r+^f+E {-T 6 m2 ^ ( e2 + iJ f - i ) v + ° (i_4) } 

where M stands for the mass of the large hole and (m, E) are the mass and orbital 
energy of the particle respectively. Note also that the force will have ^H] only negative, 
even powers of L what makes the sum quickly convergent and provides an excellent 
approximation if we evaluate the first few lower multipoles numerically. 

After this successful computation we proceeded to compute the astrophysically 
more appealing scenario of a particle in circular orbit around a nonrotating black hole. 
In this case we encountered a fundamental difficulty: The Regge- Wheeler gauge does 
not lead to C° metric coefficients at the particle location as for headon orbits. Moreover, 
in general, a Dirac's delta behaviour appear in the metric coefficients at x p . Besides, 
the regularization parameters have been computed in the harmonic gauge. 
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We have delayed until now the discussion of the important issue of the gauge 
problem. In fact, the 'MiSaTaQuWa' force (jSJ) is computed in the harmonic gauge. 
Since the 'force' is not a gauge invariant concept we need to deal with the problem of 
gauge transformations In particular, it is difficult to compute the transformation 
between a singular gauge at the location of the particle, as the Regge- Wheeler gauge, 
and a regular and locally isotropic one as the harmonic gauge. 

4. The solution to the problems of the solution 

The reason we used the Regge- Wheeler gauge was its convenience for solving General 
Relativity's field equations, i.e. reducing them to solve the Zerilli and Regge- Wheeler 
wave equations. The natural way of proceeding for a computation oriented to the 
evaluation of the self-force is to solve the ten General Relativity field equations 

+ ' ph^v = S a p, (12) 

where S a p are the source terms due to the orbiting particle, directly in the harmonic 
gauge 

g^Kp.^ = 0, where h a(3 = h a/3 - ^g a ph. (13) 

The main inconvenience to overcome are the potential numerical instabilities of this 
formulation. 

A proper use of the coordinate conditions IjlHj) that eliminates all the first time 
derivatives leading to instabilities in equations (fT2j) is possible (Barack and Lousto, in 
preparation). We can then use the algorithm of Ref. ^H] to integrate the wave equations 
with the Dirac's delta source terms appearing in f|12j) . 

The advantages of this procedure are apparent: No gauge problem, the computed 
force will be in the harmonic gauge; so are the regularization parameters. We have 
to deal with C° quantities (metric coefficients) as opposed to discontinuous waveforms 
or divergent metric coefficients, and will only need first order derivatives of them to 
compute the self-force. The procedure is in principle generalizable to Kerr backgrounds 
(although the actual implementation will not be a straightforward one). Finally, the 
non-radiative modes £ = 0,1 can be solved in practically an analytic way in the harmonic 
gauge [T3] . 

5. more problems seeking for solutions 

As we mentioned earlier the self- force is not a gauge invariant concept. A way to recover 
a physical interpretation of the corrected trajectory of the particle is to use it to compute 
the corrections to the energy and momentum radiated at infinity from the binary system. 
In order to do that we have to deal with second order perturbations of the large black 
hole background. Including this metric perturbations is needed to consistently take 
into account all 0(m/M) 2 terms (the other coming from the corrected trajectory of 
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the particle.) To second perturbative order, products of first order perturbations act as 
sources [12] . as a result, we could have singular products at the location of the particle 
(we only probed regularization of the force, not of each component of the first order 
metric and its derivatives). 

Rosenthal and Ori have recently worked out a scheme to extract the leading singular 
piece of the (quadratic) source term using the retarded field and a choice of the local 
gauge to regularize the problem. This promising approach may solve the last standing 
theoretical problem, but we know by now that new, sometimes unexpected, problems 
may still arise in our way to provide the complete solution. 
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